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TECHNICAL NOTE 3103 

COOLING REQUIREMENTS FOR STABILITY OF LAMINAR BOUNDARY LAYER 
WITH SMALL PRESSURE GRADIENT AT SUPERSONIC SPEEDS 
By George M. Low 


SUMMARY 

The amount of cooling required to stabilize the two-dimensional 
supersonic laminar boundary layer for all Reynolds numbers is calculated 
for flows with pressure gradients of a magnitude usu ally encountered 
over slender aerodynamic shapes . Only two-dimensional disturbances are 
treated in the stability calculations . 

It is determined that small pressure gradients have an appreciable 
effect on stability. The cooling due to radiation alone may suffice, at 
moderate supersonic Mach numbers, to completely stabilize the boundary 
layer over wings with favorable pressure gradients. For flows with 
adverse pressure gradients, the cooling required for complete stability 
is considerably greater than that for flat-plate flows. 


INTRODUCTION 

Because an appreciable portion of the total drag of an airplane or 
missile can be attributed to friction, a sizable reduction in friction 
drag is desirable. Such a reduction can be realized by delaying the 
transition from laminar to turbulent flow. 

Two possible explanations for the transition from laminar to turbu- 
lent flow have been advanced in the literature. Taylor (ref. l) proposed 
that the local pressure gradients accompanying a disturbance in the flow 
cause intermittent separation of the 1 amlnar boundary layer. Eddies 
formed in the separated region soon diffuse and lead to turbulent flow. 
Transition by this mechanism might be expected in high turbulence level 
wind tunnel tests or for flows where the surface roughness is appreciable. 

To llmi en and Schlichting (refs. 2 and 3), on the other hand, suggest 
that infinitesimal wave-like disturbances in the 1 amlnar layer are 
directly responsible for transition. If conditions are such that these 
infinitesimal disturbances are amplified, they will undergo increas ing 
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amounts of distortion and eventually lead to transition from 1 ami na-r to 
turbulent motion. Schubauer and Skramstad (ref. 4) have experimentally 
verified the Tollmien-Schlichting hypothesis for low speed flows . It is 
"believed that transition is caused hy laminar instability of this kind 
whenever the free-stream turbulence level is very low, and no extraneous 
disturbances such as excessive roughness exist. In free flight the tur- 
bulence level of the air is generally low, so that transition as a con- 
sequence of laminar instability would be expected whenever the aircraft 
surfaces are smooth. 

Lees and Lin (refs. 5 and 6) have developed a theory for the stabil- 
ity of the compressible laminar boundary layer based on the Tollmien- 
Schlichting concept. Their theory predicts that withdrawal of heat from 
the boundary layer has a stabilizing effect, and suggests that sufficient 
cooling will stabilize the boundary layer regardless of Reynolds number. 
Accurate and detailed calculations by Van Driest (ref. 7) based on the 
theory of Lees and Lin show that the laminar boundary layer on a flat 
plate can be completely stabilized at Mach numbers between 1 and 9. 

These results have been qualitatively substantiated by Sterpberg (ref. 

8) , who observed laminar boundary layers at Reynolds numbers as high as 
50*10 6 . 

The theory of Lees and Lin is limited to the flow over a flat plate 
(zero pressure gradient). It was subsequently shown by Laurmann (ref. 

9) and Cheng (ref. 10 ) that the criteria derived by Lees and Lin apply 
also to flows over curved surfaces provided the local velocity and tem- 
perature profiles are considered. 

In the present report calculations based on the theory of Lees and 
Lin are made for flows with small constant pressure gradients. In par- 
ticular, the cooling req uir ed to completely stabilize the laminar bound- 
ary layer for flows with pressure gradients of a magnitude usually 
encountered over thin aerodynamic shapes at supersonic speeds is calcu- 
lated. The pressure gradients are considerably smaller than the pressure 
gradients required for 1 ami nar separation. Velocity and temperature 
profiles for this calculation are obtained from reference 11. The cal- 
culations were made at the MCA Lewis laboratory during the summer of 
1953. 


ASSUMPTIONS AND LIMITATIONS 

The present work employs the results of reference 11 and therefore 
conta ins the same assumptions, discussed in detail therein, concerning 
the boundary -layer flow. These include constant specific heat, constant 
Prandtl number (Pr = 0.72), constant wall temperature, and viscosity 
proportional to temperature. Van Driest 's calculations for the stabil- 
ity of flat-plate boundary layers show that these as simp t ions do not 
lead to large errors at Mach numbers less than 3. 
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The velocity at the outer edge of the boundary layer Ug is postu- 
lated to differ only slightly from a reference velocity u r : 


lie 

— = 1 + ex 
^r 


( 1 ) 


(All symbols used in this report are defined in appendix A.) The quan- 
tity ex is small compared with unity so that the square of this quan- 
tity may be neglected. Equation (l) represents velocities which are 
increasing or decreasing linearly along the body and requires that the 
pressure gradient dp/dx be constant and equal to -reM^. Equation (l) 

x ^e 

leads to a pressure gradient parameter — = ex, which was used in 

the stability calculations . This parameter is a function not only of 
the pressure gradient, but also of the distance over which the pressure 
gradient acts. 


It should be . oted that the pressure gradients used in this report 
are the simplest case of the more general pressure distributions treated 
in reference 11. Although the stability analysis could be carried out 
for the more general pressure gradients, it was believed that the rela- 
tive effects of pressure gradients on laminar boundary layer stability 
could adequately be demonstrated for the special case of a constant 
pressure gradient. 


Because the calculations of Van Driest indicate that the assumption 
of a linear viscosity-temperature relation leads to appreciable errors 
in the stability calculations for Mach numbers greater than 3, and 
because the method of reference 11 becomes questionable at M^. > 3 for 
x due 

the values of — considered herein, all calculations were limited 

Up < 3x 

to values of M^, < 3. 

The assumptions utilized in the solution of the stability differ- 
ential equations are dis cussed in references 5 and 6. Although these 
assumptions originated from a study of the flow over flat plates, it 
was shown in references 9 and 10 that the equations derived by Lees and 
Lin apply eq ually well to flows with pressure gradients provided that 
local velocity and temperature profiles are considered in the analysis. 


STABILITY ANALYSIS 

Lees (ref. 6) suggested the possibility that, if sufficient heat is 
withdrawn from a supersonic laminar boundary layer, the flow will be 
stable at al 1 Reynolds numbers . This means that no self -excited disturb- 
ances can exist that satisfy the stability differential equations, their 
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"boundary conditions, and the physical requirement that the dimensionless 
phase velocity of the di sturbance c he greater than 1 - 

M e 

The stability differential equations are obtained by linearizing the 
equations of motion and energy for a flow with time -dependent fluctuating 
components superposed on the steady mean flow quantities. If the effect 
of viscosity is neglected in the solution of these equations, the stabil- 
ity of the fluid is found to be governed by the distribution of the pro- 
duct of density and vorticity through the boundary layer. This "inviscid" 
solution applies only at infinite Reynolds numbers. It was shown in ref- 
erence 6, however, that the effect of viscosity is destabilizing; a flow 
which is stable at Re = • may not be stable at all other Reynolds num- 
bers. In the so-called viscous solution of the stability differential 
equations, the effect of viscosity is considered to the first order; 
these solutions therefore apply at large, but finite, Reynolds numbers. 


Because the purpose of the present report is to find a wall temper- 
ature ratio below which the flow will be stable at all flight Reynolds 
numbers, the viscous solution of the stability equations is required. 
The problem reduces to finding a wall temperature ratio t w which, for 

a given free-stream Mach number and pressure gradient, simultaneously 
satisfies the following equations: 




0.580(1 + X) 
0.960X + 0.570X 2 


( 2 ) 


( 3 ) 


where 


X = 


u*'(0) It*(Tl c ) 


- 1 


u eV 


( 4 ) 


c = 1 



( 5 ) 


Equation (2) is equivalent to equation (24) of reference 6, while 
equation (3) is obtained by combining equations (l9) and (20) of ref- 
erence 6, as suggested by Bloom (ref. 12). Equations (2), (3), and (4) 
are valid in the transformed coordinate system of reference 11. A table 
of relations between the terms appearing in the present report and those 
of reference 6 is presented in appendix B. 
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Once the free-stream Mach number and pressure gradient are specified 
and a wall temperature ratio is assigned, all terms in equations (2) and 
(3) can he found from reference 11. (The relations "between these terms 
and the functions tabulated in ref . 11 are presented in appendix C . ) The 
correct wall temperature ratio satisfying both equations was found by 

trial and error.' 1 ' 


EFFECT OF PRESSURE (SiADIENT ON COOLING REQUIREMENTS FOR STABILITY 


Equations (2) and (3) were 


solved for Mach numbers between 1 and 3, 


and for values of the pressure gradient parameter imr of 


x 

u r dx 


0, +0.05, 

and +0.10. The limiting wall temperature ratio for complete stability 
t w /t is presented in figure 1 as a function of local stream Mach num- 


ber M e - At a given pressure gradient the laminar boundary layer is 
stable for al 1 Reynolds numbers if the wall temperature ratio is less 
than the value given on the curve for that particular gradient. If 
t w /t e is greater than that value, self-excited disturbances will exist 

if the Reynolds number is sufficiently high. Also shown in figure 1 are 
the wall temperature ratio for zero heat transfer and the wall tempera- 
ture ratio 5 feet aft of the leading edge of a flat plate for black-body 
radiation at an altitude of 50,000 feet. (Both the zero heat transfer 
curve and the radiation curve are altered only slightly by the pressure 
gradient.) The wall temperature ratio for radiation is obtained from a 
balance of the heat lost by the surface through radiation and the heat 
gained by the surface through convection and conduction. It was assumed, 
as is conventional, that the surface radiates to a mean receptor temper- 
ature equal to the ambient temperature. As a comparison, this heat bal- 
ance was also made for a point 2^ feet aft of the lead in g edge, and it 


was found that the resulting equilibrium wall tenperature ratio at 

M e = 3 was A percent higher than at the 5-foot station. This difference 

in temperature ratios decreases as the Mach number decreases. 


d Br. C. C. Lin and Dr. D. W. Dunn have Informed the author that 
they have obtained Improved viscous solutions of the stability differ- 
ential equations, but that calculations based on these solutions agree 
to within 2 percent with the present results for the case of zero pres- 

x du e 

sure gradient. For the case of — = 0.10 and Mg = 3.06, the new 

solutions yielded a wall temperature ratio 0.5 percent higher tha n re- 
ported herein. It is therefore believed that all the present results 
are in error by no more than about 2 percent. 
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The pressure gradients represented in this figure are of a magnitude 
that might he encountered over thin wings at supersonic speeds. For 

x du e 

example, the value of — — — at the midchord station of a 5 percent 
thick circular arc airfoil at a Mach number of 2 is 0.06. The adverse 

might be found on compression 

ramps of supersonic engine inlets . It should be noted that the parameter 
x du e 

— ^ will not be a constant along a surface when the pressure gradient 

r 

is constant. Thus a given stability limit, an plotted in figure 1, 
applies only at one point on a wing, and the boundary layer may be more 
stable (or less stable) at other chordwise stations. 


/x du e \ 

pressure gradients I — — — negative^ 


stream Mach number of 2, and for a value of 


equal to 0.10, the 


It is evident that the effect of a small pressure gradient on the 
cooling requirements for complete stability is appreciable. At a local 

x_ ^e 

u r dyi 

wall temperature ratio for stability is not far below the zero-heat- 
transfer temperature ratio. For this pressure gradient the cooling due 
to radiation alone is seen to stabilize the boundary layer for Mg> 1.7. 
For flows with adverse pressure gradients, the wall temperature ratio 
for complete stability is considerably lower than for flat-plate flows. 


The effects of pressure gradient may be demonstrated in terms of the 
relative rate of heat transfer required, in the absence of radiation, to 
stabilize flows with pressure gradients as compared with flat-plate flows. 
The ratio of heat-transfer rates q/q^p, as obtained from the temperature 
ratios of figure 1 and from equation (63) of reference 11, is plotted as 
a function of local stream Mach number in figure 2, and applies at an 
ambient air temperature of -67° F. At a Mach number of 2.5, only 5 per- 
cent of the cooling required to stabilize the flow over a flat plate 

x du e 

will suffice to stabilize the flow over a wing when — — ^ = 0.10. 

For the most adverse pressure gradient considered, on the other hand, 
the ratio q/q^p reaches a peak value of about 2.2. 

There appears to be an optimum Mach number, in the neighborhood of 2, 
where pressure gradients have a large effect on the cooling requirements 
for complete stability. At higher Mach numbers the frictional heating 
terms in the bo undar y layer energy equation became of greater importance, 
and the relative effect of pressure gradient becomes less important. 
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A possible shortcoming of the present calculations, in common with 
the shortcoming of previous papers, is that only two-dimensional disturb 
ances are treated. It was shown by Squire (ref. 13) that, for incompres 
sible flow, two-dimensional disturbances are always more destabilizing 
than three-dimensional disturbances . Since the present calculations 
have been completed, however, it has been reported by Dunn and Lin 
(ref. 14) that Squire's theorem cannot be extended to supersonic flows, 
and that under certain conditions three-dimensional disturbances will 
be more destabilizing than two-dimensional disturbances. The trends of 
the two-dimensional theory have, however, been verified experimentally 
(ref. 8). It is believed, therefore, that the present results at least 
qualitatively describe the effects of pressure gradients on boundary 
layer stability. 


CONCLUDING REMARKS 


The amount of cooling required to completely stabilize a two- 
dimensional laminar boundary layer with small pressure gradients in 
Supersonic flow has been calculated. It was found that the wall temper- 
ature ratio for complete stability approaches the zero-heat -transfer 
temperature ratio for reasonably small favorable pressure gradients. 

At a local stream Mach number of 2.5, for example, only 5 percent of 
the cooling required to stabilize the flow over a flat plate will suf- 
fice to stabilize the flow over a wing when the pressure gradient param- 

jj. dUg 

eter ^ equals 0.10. Under these conditions the cooling due to 

radiation may be adequate to stabilize the boundary layer for all 
Reynolds numbers. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 3, 1953 
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APPENDIX A 
SYMBOIS 

The following symbols are used in this report: 
c phase velocity of disturbance divided by u g 
M Mach number 

q local rate of heat transfer 

t static temperature 

t* t/t r 

u velocity in x-direction 

u* u/u^ 

V function appearing in eqs. (2) and (3) 
x distance along surface measured from leading edge 
X ratio of specific heats 

e RTTioii quantity - measure of shape and magnitude of velocity dis- 
tribution at outer edge of boundary layer 

T) characteristic variable (see appendix B) 

\ wave length of disturbance 

Subscripts: 

u 

c value of function when — t=- c 

u e 

e conditions at outer edge of bou n da r y layer 
FP equivalent flat -plate value 

r reference condition (see ref. 11 ) 


w 


c ondi ti ons at wall or surface 
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Superscripts : 

' differentiation with respect to rj 

* dim ensionless quantity 


Special notation: 


The symbol I preceding a quantity indicates integration from zero 
to r\j for example. 





t*C?)d£ 


CM 

I 
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APPENDIX B 


RELATION BETWEEN NOTATION OF REFERENCE 6 AND NOTATION 
EMPLOYED IN PRESENT REPORT 

Because the present stability calculations are based on the velocity 
and temperature distributions of reference 11, it was necessary to obtain 
the stability equations in terms of the variables used in that reference. 
A table of equivalent relations which are required in order to obtain 
equations (2), (3), and (i) from the equations of reference 6 is there- 
fore presented: 


Lees' notation 
(ref. 6) 

Present notation 

c 

c 

T 

t»/tf 

V 

V 

w 

u^/u* 
1 e 

y 

It*(n) 

X 

X 

, a 

, d d _ 1 d 

“ ^y 

d V ay t* dT i 

„ i 2 

.. _ d 2 d 2 _ 1 d 2 t* f d 

ay 2 

dq2 ay2 (t*)2 dT|2 (t*)^ ^ 

Subscript 1 

Subscript w 

Subscript 0 

Subscript e 
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APEEHDIX C 


DEFINITION OF TERMS USED IN STABILITY CAICULATION 

The terms required in the stability calculation for flows with pres- 
sure gradients are related to the functions tabulated in reference 11 
as follows: 

\ f (n) + ex g*(ij) 

U*’ = \ f"(Tl) + 6X g"(T}) 

u*" = \ f-Cn) + ex S 1,1 (tj ) 
t* » 1 + K s(ti) + ^ Mj. 2 t( n ) - 2 ex H(tj)| 
t*’ = K s'(ri) + [r.(,) - 2ex H’Ol)] 

It*(ri) = T] + K Is( n ) |lr(ti) - 2ex IH(ij)j 


where 


s(n) = SnOl) + Mr 2 gi2(n) + KBijCtj) 

H(t] ) - SuCt) + ^r 2 H 12 ("H ) + ^13 (t) + H i 40 l) + H 15 (ii) 

Mj. Mj,* 


and 


K “ iW ' 1 " ^ “r 2 


All functions of t] with the exception of f "'(tj) and g"'(T]) 
are tabulated in reference 11. The function g'"^) is presented in 
table I of the present report, while f is obtained from the Blasius 
equation 

f'»= - ff" 
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The relation between the local stream Mach number Mg and the refer- 
ence Mach number M,, is obtained from equation (l) and the corresponding 
equation for the local stream temperature 

t* = 1 - (r-l) Mp 2 sx 

Thus 


M e =M r 


1 + ex 
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TABLE I. - THE FUNCTION OF g|” (t|) 




a t »» 
*12 


0 

- 4.0000 

- 0.6781 

- 8.2990 

.1 

- 4.0000 

-.6694 

- 7.8084 

.2 

- 3.9996 

-.6433 

- 7.3179 

.3 

- 3.9972 

-.6001 

- 6.8275 

.4 

- 3.9887 

-.5404 

- 6.3352 

.5 

- 3.9670 

-.4656 

- 5.8374 

.6 

- 3.9219 

-.3779 

- 5.3284 

.7 

- 3.8409 

-.2807 

- 4.8016 

.8 

- 3.7105 

-.1781 

- 4.2507 

.9 

- 3 .5178 

-.0755 

- 3.6715 

1.0 

- 3.2538 

.0214 

- 3.0637 

1.1 

- 2.9147 

.1068 

- 2.4326 

1.2 

- 2.5047 

.1757 

- 1.7898 

1.3 

- 2.0358 

.2241 

- 1.1526 

1.4 

- 1.5283 

.2501 

-.5428 

1.5 

- 1.0080 

.2539 

.0163 

1.6 

-.5037 

.2378 

.5023 

1.7 

-.0432 

.2060 

.8973 

1.8 

.3501 

.1636 

1.1902 

1.9 

.6597 

.1163 

1.3781 

2.0 

.8780 

.0693 

1.4661 

2.1 

1.0060 

.0270 

1.4663 

2.2 

1.0522 

-.0078 

1.3957 

2.3 

1.0306 

-.0336 

1.2738 

2.4 

.9583 

-.0503 

1.1202 

2.5 

.8527 

-.0587 

.9525 

2.6 

.7299 

-.0603 

.7850 

2.7 

.6033 

-.0571 

.6284 

2.8 

.4828 

-.0508 

.4894 

2.9 

.3748 

-.0430 

.3712 

3.0 

.2828 

-.0348 

.2746 

3.1 

.2074 

-.0271 

.1981 

3.2 

.1485 

-.0204 

.1400 

3.3 

.1028 

-.0147 

.0960 

3.4 

.0697 

-.0103 

.0647 

3.5 

.0463 

-.0071 

.0427 

3.6 

.0298 

-.0047 

.0276 

3.7 

.0189 

-.0030 

.0175 

3.8 

.0116 

-.0020 

.0108 

3.9 

.0071 

-.0012 

.0067 

4.0 

.0041 

-.0006 

.0040 

4.1 

.0024 

-.0004 

.0024 

4.2 

.0014 

-.0002 

.0014 

4.3 

.0007 

-.0001 

.0008 

4.4 

.0004 

-.0001 

.0004 

4.5 

.0002 

0 

.0002 




Kail temperature ratio, tv/W 


Oil/ 



K 


Figure 1. - Limiting wall toaperature required, for oomplete etatllizatioa of lam -I nay boundary layer. 
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Relative heat-transfer rates, q/qjrp 
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1.0 1.4 1.8 2.2 2.6 3.0 3.4 

Local stream Mach number. Mg 


Figure 2. - Cooling requirements for flows with pressure gradients 
referred, to cooling requirements for flat-plate flows. 
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